Quantized ionic conductance in nanopores 
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Ionic transport in nanopores is a fundamentally and technologically important problem in view of 
its occurrence in biological processes and its impact on novel DNA sequencing applications. Using 
microscopic calculations, here we show that ion transport may exhibit strong nonlinearities as a 
function of the pore radius reminiscent of the conductance quantization steps as a function of the 
transverse cross section of quantum point contacts. In the present case, however, conductance steps 
originate from the break up of the hydration layers that form around ions in aqueous solution. 
Once in the pore, the water molecules form wavelike structures due to multiple scattering at the 
surface of the pore walls and interference with the radial waves around the ion. We discuss these 
effects as well as the conditions under which the step-like features in the ionic conductance should 
be experimentally observable. 



Over the last decade there have been tremendous ad- 
vances in both the fabrication of nanopores and their 
use for molecular recognition and nucleic acid analy- 
sis [ij. Experimental characterization of molecules has 
mostly relied on measuring changes in the ionic current 
through the pore [1, H, 0, 0| , but other ways to probe 
single molecules in nanopores may come from embed- 
ding nanoscale sensors within the pore or nanochannels 
0, i, i, E, [HI, [H, El. However, electronic fluctu- 
ations due to the dynamical ionic and aqueous environ- 
ment will affect the type of signals and noise these sen- 
sors detect. Therefore, understanding the electrostatics 
of ions in water at atomic length scales is crucial in our 
ability to design functional single-molecule sensors and 
to interpret their output, and will also provide new in- 
sight into the operation of the ubiquitous biological ion 
channels. 

Several studies have examined the electrostatics of ions 
in nanopores using continuum models for the dielectric 
properties of water 14, E3, E3, 12 1 . Within a contin- 



uum model, the nanopore electrostatic environment is 
essentially one-dimensional due to the large difference 
of dielectric constants between water and the surround- 
ing pore material [3]. Thus, according to these models, 
there is a large electrostatic energy penalty to move an 
ion from the exterior of the pore to its interior 14 , 11 , 
with small amounts of surface charge able to drastically 
reduce this energy penalty [31 . Although continuum 
models can highlight some generic features of ionic cur- 
rents and blockades, such as the effect of surface charges, 
they miss important effects related to the microscopic 
physical structure of water molecules around ions. 

Here, we examine ionic transport from the point of 
view of these nanoscale features (see schematic in Fig. 
(HI)). Using all-atom molecular dynamics simulations, we 
calculate the structural and electrostatic details of a sin- 
gle ion in an aqueous environment both in and out of a 
cylindrical nanopore (see Methods for details). Ions in 
solution create local structures in the surrounding wa- 



Figure 1: Schematic of a single anion with its first hydration 
layer that moves into a nanopore. Due to spatial constraints 
many water molecules need to be stripped off for the ion to 
enter the pore, creating nonlinearities in ionic conductance. 



ter, known as hydration layers [18|, and characterized by 
oscillations in the water density. The length scales asso- 
ciated with these oscillations and their decay are in the 
range 2 — 10 A, i.e., comparable to the radial dimensions 
of some nanopores. In addition, we show that when the 
ion is inside the pore, water molecules form wavelike fea- 
tures due to both the interaction with the walls and the 
water molecules in the partially broken hydration lay- 
ers. Starting with these microscopic considerations, we 
propose a model system based on the concept of energy 
shells, which incorporates the hydration layer structure 
into the energetic barrier created by the pore. We pre- 
dict that the ionic conductance of a nanopore manifests 
step-like features as a function of the pore radius similar 
to the quantized conductance of quantum point contacts 
as a function of their effective cross section (see, e.g., 
[l^). Finally, we examine the influence of noise on ionic 
transport and discuss the experimental conditions under 
which the above step-like features are most likely to be 
observed. 

To be specific, we focus our attention on an isolated 
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chlorine ion but similar considerations will apply to other 
types of ions and situations. Immediately around the 
CI anion, the water molecules orient themselves so that 
a single hydrogen from each molecule points inward to- 
ward the ion (see Fig. (jT])). The distribution of water 
molecules forms into layers as shown in Fig. The 
inner most layer is very tightly bound, and subsequent 
layers are spaced at 2.0 — 2.3 A. These findings are in 
good agreement with neutron diffraction and X-ray ab- 
sorption experiments that measure the radial distribution 
of water [20|, which give a Cl-0 peak at '--^ 3.1 A, as well 
as molecular dynamics simulations performed with dif- 
ferent force field parameters (2]] |. In addition, we can 
acquire further information not directly accessible from 
experiments, such as the microscopic electric field due to 
the ions and water (shown in the inset of Fig. (Il])(b)). 
This field shows oscillations corresponding to the hydra- 
tion layers and is similar to a series of alternating charged 
surfaces. A continuum picture, however, would not cap- 
ture these microscopic details that are responsible for the 
conductance steps we predict. 

If the anion is now placed in the nanopore, the hydra- 
tion layers are affected. For large pores, the hydration 
layer structure is identical to that in bulk water. How- 
ever, as the pore is taken to nanoscale dimensions, even- 
tually the pore walls (or effective pore wall, which may 
be, e.g., a layer of tightly bound water molecules) will 
force the hydration layer to be partially broken because 
it can not fit within the pore. Figures [S^c, d) show this 
effect. Within a 15 A radius pore, the first and second 
hydration layers are still present. On the other hand, 
for a 12 A pore radius, the first layer is still intact, but 
the second layer is partially broken. Indeed, it is almost 
a spherical shell cutoff by the effective pore wall which 
comes at about 5 A from the axis of the pore. In order 
to gain physical insight, we propose the following model 
system that captures this essential aspect of the presence 
of hydration layers around the ion. 




Figure 2: Water density around a CI anion both in bulk and 
inside cylindrical nanopores of different radii, (a) Water den- 
sity (using the oxygen as the location of the water molecule) 
around a CI anion in bulk water, (b) Water density versus ra- 
dial distance from a CI anion. At distances greater than 10 A 
the water density is approximately 0.036 molecules/ A"^, which 
is slightly more than bulk water density of 0.033 molecules/A^ 
at room temperature. Some edge effects are present around 
15 A. The distribution of water molecules shows oscillations 
corresponding to the formation of hydration layers, which is 
also reflected in the electric field. The inset shows the time 
averaged radial electric field (versus the radial distance) from 
both the ion and water dipoles (red line) and from just the 
ion (black dashed line), (c) Water density around a CI anion 
inside a 15 A radius pore (the pore walls are at the edge of 
the figure). There are wavelike features that are due to inter- 
ference patterns between oscillations reflecting off the walls of 
the pore (that set an effective pore radius) and those around 
the ion. However, the first two layers around the anion are stiU 
present, (d) Water density around a CI anion inside a 12 A 
radius pore (the pore walls are indicated by vertical, thick 
black lines). Here the effective radius is approximately 5 A, 
thus the second hydration layer interacts with the effective 
pore wall and it is partly broken. 



Model System - The hydration layers stay partially in- 
tact except for the portion of the hydration surface cut off 
by the pore walls, i.e., for layers with radii greater than 
the effective pore radius. In other words, only the por- 
tion of the layer along the axis of the pore remains intact. 
We consider then a set of surfaces placed at each hydra- 
tion layer, z, at radius Ri. Each surface represents the 
area where the water dipoles fluctuate, givi ng the time- 
averaged dipole layers. Ignaczak et al. [23] have found 
that for small water clusters around a CI anion, the inter- 
nal energy of each water is approximately linear in the 
number of water molecules in the cluster. Since inter- 
actions (van der Waals and electrostatic) with a low di- 
electric pore are small compared to water-ion and water- 
water interactions , this suggests writing the internal 
energy contained in a partially intact hydration layer as 



u^ - m (1) 

where fi is the fraction of the layer present in the pore 
and U° is the energy difference of that stored in the intact 
water layer and the water in bulk. 

In our model, the energy barrier is due to the stripping 
off of a fraction of the layer, fi, i.e., the fraction of a 
spherical surface at Ri that remains within the pore of 
effective radius Rp. When the ion translocates along the 
pore axis, the surface area that remains in layer i is given 

by 

Si = 29 {R^ - Rp) / d(p dORj sixvO (2) 
Jo Jo 

where Q (x) is the step function and 6c — sin^^ Rp/ Ri. 
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When R„ < Ri, the fraction of the surface left is 



fi (Rp) 



(3) 



The internal energy barrier as function of pore radius is 
then given by 



A;7(i?p) = ^(/.(i?p)-i)[/°, 



(4) 



where the summation is over the layers. The free en- 
ergy change of hydration is dominated by the contribu- 
tion from the total internal energy contained within the 
hydration layers, J^i U", plus a bulk dielectric contribu- 
tion from water outside the layers. Previous calculations 
give the total internal energy change in the range -3.5 to 
-3.7 eV fii, '24', '25^. The form of the electric field shown 
in Fig. m^b) suggests that the time-averaged microscopic 
distribution of waters can be viewed as a set of spheri- 
cal layers of alternating charge. Thus, the energy stored 
within each layer should be well estimated by using a 
Born solvation calculation, which, by construction, will 
sum to the total internal energy contribution . The phys- 
ical content of this calculation is that of the solvation of 
embedded "quasiparticles" : the total solvation energy is 
given by solvating the ion; the energy of the first layer 
is the difference of solvating the ion and solvating the 
ion plus the first layer, i.e., solvating a "quasiparticle" 
consisting of the ion and a cluster of waters. With this 
picture in mind, the energy within each layer is given by 
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where k is the dielectric constant of water and i?| are 
the inner (outer) radii demarcating the hydration layer, 
and we have used the effective ion charge of le, which is 
consistent with the field recovering its bare ion value after 
each hydration layer (see inset of Fig. [H^b)). From the 
water density oscillations shown in Fig. [21 these radii are 
2.0 and 3.9 for the first layer; 3.9 and 6.2 for the second; 
6.2 and 8.5 for the third. The total solvation energy is 
then -3.6 eV and the layer energies are -1.7, -0.7, and - 
0.3 eV for the first, second, and third layers, respectively. 
Immediately, we can understand why the third layer is 
absent within the pore: van der Waals interactions 23 1 
and entropic contributions from the water structure 26 1 
are of the same magnitude. 

In addition to the internal energy change, there is also 
an entropic contribution to the free energy change that 
comes from removing a single ion from solution and lo- 
calizing it in the pore. Assuming an ideal ionic solution, 
this entropic contribution is AS" = fcs In (Vpn), where Vp 
is the volume of the pore, n is the bulk ionic concentra- 
tion, and ks is the Boltzmann constant. The final free 
energy difference, AF = AU — TAS, is plotted in Fig. 




Figure 3: (a) Free energy barrier to bring an ion into a 
nanopore and the ionic current at room temperature as a 
function of the effective pore radius. The dotted red line in- 
dicates the ionic current without a free energy barrier and 
the dashed blue line is the current with just the entropic 
barrier. The current is calculated at 120 mV and 1 M 
bulk salt concentration assuming the ionic current is given 
by / = 2TvRpefinE, where the density n is proportional to 
a Boltzmann factor containing the energy barrier, namely 
n — no exp{—AF/kT), with no the bulk salt concentration, 
fi the ion mobility, and E the strength of the electric field. 
A typical ion mobility is /i = 2.28 X IQ-^m^/Vs (see, e.g., 
Refs. and p^). (b) dl /dRp as a function of pore radius 
and noise strength. For zero noise, the step-like features in 
the ionic current give rise to well defined peaks in the deriva- 
tive of the current. However, as the noise is increased, the 
peaks smooth out. This occurs beyond a noise strength of 
5 = 0.2 — 0.3 A for the second and third layers (see text for 

details), (c) Relative noise, Alrei = {{l^)-{lfY^'^/{I), 
versus pore radius for the outermost hydration layer and for 
C = 0.05, 0.15, 0.25, 0.35, 0.45. If the pore radius is tuned 
to the radius of a hydration layer, fluctuations give rise to a 
peak in the relative noise, which gives an alternative method 
to observe the effect of the hydration layer. 



Q . This is the main result of this work: hydration layers 
hold energy in a shell-like structure that causes step-like 
features in the energy barrier to ionic transport. These 
step-like features generate corresponding features in the 
ionic current shown in Fig. ([3]). 

For the remaining portion of this paper, we discuss the 
conditions under which these features can be observed ex- 
perimentally and address how the assumptions we have 
made affect our conclusions. First, we have assumed that 
only the single ion barrier is important, i.e., that there 
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are no correlated transport processes. We can obtain 
an indication of whether experiments are consistent with 
this assumption by considering the ionic concentration. 
If the bulk salt concentration is 1 M, a typical ion (count- 
ing both anions and cations) will be contained in a vol- 
ume with linear dimension 9.4 A. This is large enough 
to almost house both the first and second hydration lay- 
ers. Since in the absence of significant surface charge 
on the pore walls, one expects the ionic concentration 
within the pore to be lower, a bulk concentration of 1 
M is consistent with having ions translocate through the 
pore independently. AIM concentration is likely to be 
small enough to allow at least the effects of the first ionic 
layer to be observed, i.e., a sudden drop in the current 
when the pore radius is comparable to the first hydration 
layer. For arbitrary concentrations, the ion-ion distance 
goes as ~ 9.4/no^^A. If the ion is to be contained in a 
linear dimension of 20 A to allow all hydration layers to 
be present, then one would want a concentration of about 
0.1 M. Low concentrations are therefore key to observe 
the outermost hydration layers. Furthermore, the ionic 
concentration itself may be used to observe the step-like 
features. Namely, at certain concentrations the anion- 
anion (cation-cation) distance would be small enough so 
that the hydration layers should break with increasing 
concentration, and a corresponding nonlincarity would 
be observed. 

Second, we have neglected energetic fluctuations due to 
other effects present, such as thermal noise, ion-ion inter- 
action, rough pore surfaces, a distribution of ion paths 
through the pore, and the fact that hydration layers are 
not exactly represented by their time-averaged oscilla- 
tions (therefore, some ions with intact hydration layers 
could go into pores smaller than the radius Ri would 
predict). In order to determine how robust the step-like 
features are to these noise sources, we add in our calcu- 
lations Gaussian noise in the pore radius. In particular, 
we define a noise strength parameter ^ (i.e., the stan- 
dard deviation of the Gaussian noise) to describe both 
energetic and path fluctuations. In Fig. ([2]), we plot 
dl/dRp as a function of Rp and For no noise, ^ = A, 
there are well defined peaks in dl/dRp. However, at a 
noise strength of f = 0.2 — 0.3 A, the peaks from the sec- 
ond and third layers are smoothed out. Given the width 
of the hydration layers as seen in Fig. the effects 

from the second and third layers may not be observable 
when transport is due to ions like Cl~, or K"^. However, 
other ions, such as the divalent Mg or trivalent Al ions 
[29L [30I I , as well as some monovalent ions jl^ , have more 
strongly bound second and third layers. Together with 
using tunable properties, such as temperature [2^ Isot . 
this may allow observing the step-like features in the cur- 
rent. In addition, the noise gives an alternative way to 
observe the effect of the hydration layers. As shown in 
Fig. ©(c), the relative noise peaks when the pore ra- 
dius is near that of a hydration layer. This effect can be 



understood in terms of a two-state picture: fiuctuations 
that decrease the pore radius energetically "close" the 
pore, whereas fluctuations that increase the pore radius 
"open" the pore, essentially creating an instability in the 
current. Such an effect could also be induced artificially 
by sampling over a distribution of pore radii centered 
around a layer radius. The other features of the relative 
noise can also be understood within this picture. 

Third, we have assumed ion translocation is along the 
axis of the pore. In reality, the ions can move off center. 
There are different factors which cause this. First, the 
hydration layers at any given instant are not spherically 
symmetric. In fact, some ions have asymmetric hydration 
layers where at a given time, for instance, the first layer 
is offset to one side of the ion [2^. However, this fact 
cannot destroy the effect we predict: in these cases, the 
hydration layer may be able to "squeeze" into a smaller 
radius pore by moving off center, but eventually it has 
to break in a highly nonlinear way. Another cause, how- 
ever, is something that our model predicts: if the energy 
is contained in a spherical layer, then the ion actually 
tends to move off center. This is a consequence of geom- 
etry: more of the spherical surface can fit into the pore 
if the ion moves toward the pore wall. However, the ion 
cannot move too far from the center, otherwise some of 
the inner layers would have to be partly destroyed with 
consequent energy penalty. Neither of these two causes 
would change our predictions, although they may change 
some details, such as at what pore radius the conduc- 
tance jumps occur. Finally, a strong affinity of the ion 
for the pore wall may occur, for instance, if the pore were 
significantly charged. This effect may indeed reduce the 
energy barrier for ionic transport considerably, making 
all but the first hydration layer breakup undetectable in 
the ionic current. The charging of the pore walls, which 
depends heavily on pore material and treatment [3l|, as 
well as other factors, can again be somewhat resolved by 
observing transport with different ion and pore types. 

Before we conclude, we briefly mention the other water 
patterns within the pore. These patterns depend on pore 
size - different size pores can give maxima and minima in 
the water density along the pore axis (3^ . Further, we 
find that when the pore radius is 10 A (effective radius 
3 A) there is a transition in the water structure of the 
Carefully controlled ionic transport 
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pore 

measurements may also be able to detect these changes 
in addition to the nonlinear features due to the hydration 
layers. 

Conclusions - We have examined ionic transport 
through nanopores using both molecular dynamics and 
electrostatic calculations. A simple form of the free en- 
ergy barrier can capture an important effect of hydration 
layers - namely that the layers produce step-like features 
in the energy barrier for ions to enter a nanopore, and 
thus also the ionic conductance, as a function of pore ra- 
dius. This effect is the classical counterpart of the elec- 
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tronic quantized conductance one observes in quantum 
point contacts as a function of their cross section. Like 
the quantum case, one needs sufficient experimental con- 
trol of the nanopore characteristics in order to observe 
this effect. Irrespective, a drop in the current should 
be observable when the first hydration layer is hit, sim- 
ilar to what is thought to occur in the open and closed 
states of some biological pores [H, [s^, IS] • However, un- 
der the right conditions, multiple current drops should 
be observed. Further, the relative noise has a charac- 
teristic peak when the pore radius is near a hydration 
layer radius, and this gives an alternative way to observe 
the effect of hydration. Due to the growing importance 
of ionic transport in DNA sequencing, this effect may 
have consequences for future technologies. In addition, 
the operation of biological ion channels has many fac- 
tors associated with it, such as favorable internal sites 
[3^ that compensate for the loss of hydration (favorable 
interactions can also exist in artificial pores 4^, 4l|) or 
correlated ionic processes [42|. Thus, investigations with 
more well-controlled synthetic pores can help shed light 
on the different factors at play and we hope our work will 
motivate experiments to specifically examine the role of 
dehydration. 

Methods - The molecular dynamics simulations are 
performed with the NAMD2 packag e_j43l | and the sim- 



ulation details are as in Refs. Y, 3, l9|. The TIPS [44 1 
model is used for the description of water, which is known 
to reproduce the dielectric properties of bulk water ■ 
For all simulations, a single chlorine anion is taken as 
fixed since its thermal wavelength at room temperature is 
- O.2A. The cyhndrical, Si3N4 nanopore surfaces carry 
a charge of 6 — 8 x 10"'' e/atom for the 15 A and 12 A 
radius pores. The density plots for an ion in a spherical 
droplet of water are calculated by averaging the last 900 
out of a 1,000 snapshots (taken at 1 ps intervals) within 
spherical shells of thickness 1 A. For an ion inside a pore, 
the last 4000 out of 5000 snapshots (taken at 0.1 ps in- 
tervals) are averaged within cylindrical shells of thickness 
1 A in both the radial and vertical directions. The elec- 
tric field is calculated by averaging the ion and dipole 
moment contributions over the snapshots. 
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